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Abstract 

We construct a separation of variables for the classical n-particle Ruijsenaars system 
(the relativistic analog of the elliptic Calogero-Moser system). The separated coor- 
dinates appear as the poles of the properly normalised eigenvector ( Baker- Akhiezer 
function) of the corresponding Lax matrix. Two different normalisations of the BA 
functions are analysed. The canonicity of the separated variables is verified with 
the use of r-matrix technique. The explicit expressions for the generating function 
of the separating canonical transform are given in the simplest cases n = 2 and 
n = 3. Taking nonrelativistic limit we also construct a separation of variables for 
the elliptic Calogero-Moser system. 
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1. Introduction 



One of the most powerful methods in studies of Liouville integrable systems is that 
of Separation of Variables (SoV). Originated with the development of the Hamilto- 
nian mechanics as a method to solve the Hamilton- Jacobi equation for particular 
Hamiltonians, nowadays it has been applied to many families of finite-dimensional 
(Liouville) integrable systems (see recent review fl3lH ). 

For a very long time a great deal of attention has been given to so-called co- 
ordinate separation of variables or to separation in the configuration space (see, 
for instance, |8], |28|, [14], [|, |], |1(], |31[ and references therein). In this case the 



separation variables Uj do not depend on the momenta pi and are functions of the 
coordinates Xi only: 

Uj = Uj(xi, . . .,x N ) . 

Such kinds of integrable systems admitting a coordinate (local) separation of vari- 
ables were studied in detail, although in the same time it was understood that far 
not every Liouville integrable system can be separated through such a transition to 
new "coordinates" Uj. The class of admissible transformations should be enlarged 
for a generic integrable system up to a general canonical transformation 



u 



j 



Uj(xi, . . . ,X N ,Pi, . . . ,p N ) , Vj = Vj(xi, ...,x N ,px,.. .,Pn) ■ 



In the context of the Inverse Scattering Method ^ |3l| the separation vari- 
ables (u, v) appear usually as pairs of canonically conjugate variables sitting on the 
spectral curve of the related n x n Lax matrix L(u). The coordinates Uj are ob- 
tained respectively as the poles of the associated Baker- Akhiezer (BA) function /(it) 
satisfying the linear problem 

L(u) f(u) = v f(u) , /(«) = (fx(u), . . . , f n (u)Y , 

with some fixed normalisation a (it) 

n 
3=1 

The method of SoV in such a formulation was successfully applied to many particular 
integrable systems, here are some of the relevant references |25|, ^6], |27|, [j], ^4], ^T], 

0,0, H, MM- 



In the present paper we prove the SoV for the classical n-particle Ruijsenaars 
system with the n x n Lax matrix found in p3j and with the Hamiltonian 



j=i 



a(xj - x k ) 



where o~{x) is the Weierstrass cr-function, A G M is a parameter of the model and 
(pj, Xj) are canonical Darboux variables. It is shown that the method of SoV applies 
to this system if we use the standard normalisation vector a 

a = d = (0,0, . . .,0, 1) , i.e. f n (u) = l. 
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The structure of the paper is the following. In Section 2 we collect known 
information about the Ruijsenaars system (Lax matrix, integrals of motion, etc). In 
Section 3 we give an overview of the method of separation of variables and apply it 
then, in Section 4, to the system in question. In that key Section we also discuss the 
possibility of an alternative choice for the normalisation vector a(u). The generating 
functions of the canonical separating transform given in terms of the initial and 
separation variables are constructed in Section 5 in explicit form for the case of two 
and three degrees of freedom. We also provide the separation of variables for the 
nonrelativistic limit A — > to the elliptic Calogero-Moser system in Section 6. The 
Section 7 contains some concluding remarks. 



2. The system 

Let us first recall some properties of the Weierstrass functions which we will need 
in the main text. Let 2u\ t 2 G C be a fixed pair of the primitive periods and T = 
2o>iZ + 2a;2Z the corresponding period lattice. Let us fix also the primitive domain 
V := {z = 1ui\x + 2uj 2 y I x iV £ [0)1)} sucn that V ~ C/T. The Weierstrass 
sigma-function is defined by the infinite product (cf., for instance, 

(2.1) 



a(x) = x I] (l-f) exp^ + i(^) 2 
7 er\{o} 

the relations between a-, (- and p- functions being given by 

COO = ^ > = -C'(s) , (2-2) 

a[x) 

where o~(x) and ((x) are odd functions and p(x) is an even function of its argument. 
We recall also that the o~(x) is an entire function, and ((x) is a meromorphic function 
having simple poles at uiki, both being quasi-periodic, obeying 

C(x + 2u lj2 ) = C(x) + 2 77 1>2 , o{x + 2u; li2 ) = -o{x) e 2r >^ x+ ^ , 

in which 771,2 satisfy tji 0J2 — i]2 oj\ = y , whereas p{x) is doubly periodic. From 
an algebraic point of view, the most important property of these functions is the 
existence of a number of functional relations, the most fundamental being 

((a) + C(/3) + C(7) ~ C(o! + P + 7) = T\ TW\ ( \ ( \ 2 " 3 

a (a) a(fj) (7(7) a (a + (J + 7) 

which can be cast into the following form 

$ K (x) $ K (y) = $ K (s + y) [C(k) + C(x) + C(y) -((K + x + y)} (2.4) 
with the function $ K (x) defined as follows: 

a(x + k) 



a(x\ <7{K) 
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Two other useful identities have the form 



^n-k(a - b) + b) + a) - - y) ® K (y + a) + b) (2.5) 

= $ K (x + a) <$>~ K {y + b) [ ({a - b) + C(x + b) - ({x -y)-({y + a)}, 

$> K - ii (x-y)$ K (y + a)<$>~ K (x + a) (2.6) 
= $ K (s + a) <& h {y + a) [ ((x - y) - C(k + z + a) + £(« + y + a) + ((k - k) ] . 

The generalised Cauchy identity has the following form || 

det («.(„ -,,)) = ».(£) .(E) ° jfe ~ *' ^ - »*> (2.7) 

IIfc,i - 2//) 

where E = E»(^i - 

Now we can introduce the n-particle (Ai-i type) Ruijsenaars system p3| . It is 
an integrable system with the following integrals of motion (i — 1, . . . , n) 

*.= e «p(e») n g( rV - (28) 

JC{l,...,n} / JSJ U l A J ^k) 

\J\=i v 7 fce{i,...,n}\j 

The variables (pj,Xj), j = l,...,n, on a 2n-dimensional symplectic manifold form 
a canonical system, i.e. they possess the Poisson brackets 

{pj,x k } = {xj,x k } = , {pj,x k } = 5 jk , j,k = l,...,n, (2.9) 

or, equivalently, the symplectic form uj is expressed as u = J2j dPj A dxj = d(J2j Pj 
dxj). The A is a parameter of the model. This system was proposed by Ruijsenaars 
as a relativistic analog of the Calogero-Moser system. 

Proposition 1 ( [|23|| ) . The Hamiltonians Hj Poisson commute 

{H j7 H k } = 0, j,k = l,...,n. (2.10) 



The Lax matrix for this model has the form 

n 

L ( U ) = hi ®u(%i ~ X 3 + ^) E ij > h i '■ = 

where the matrix have the following entries: (E, 



n 



a [Xi — x 



A) 



cr[Xi 



(2.11) 



ij ) kl 



S ik Sji. Notice that Rui- 
jsenaars |23| used another gauge of the momenta such that two are connected by 
the following canonical transformation: 



p, -> Pi + log n , 

\ o-ixi-Xj - A) 



a(xi - xj + A) 



Xj 



X i 



(2.12) 
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Proposition 2 ([|23|j). The characteristic polynomial of the matrix L(u) ($.11 ) 



generates the Hamiltonians $Zl\) 



where we assume Hq = 1. 



3. The method 

Recall, first, the standard definitions of Liouville integrability and SoV in the Hamil- 
ton- Jacobi equation 0. An integrable Hamiltonian system with N degrees of free- 
dom is determined by a 2A r -dimensional symplectic manifold (phase space) and iV 
independent functions (Hamiltonians) Hj commuting with respect to the Poisson 
bracket 

{H j} H k } = 0, j,k = l,...,N. (3.1) 

To find a SoV means then to find a canonical transformation M : (x,p) i— > (u, v), 
M : Hi(x,p) i— > Hi(u, v) such that there exist N relations 

$ j (u j ,v j ;H 1 ,...,H N ) = 0, j = 1, . . . , N , (3.2) 

separating the variables Uj. The most common way to describe a canonical trans- 
formation is the one in terms of its generating function F{u\x). 

Presently, no algorithm is known of constructing a SoV for any given integrable 
system. Nevertheless, there exists a fairly effective practical recipe based on the 
classical inverse scattering method. A detailed description of the procedure with 
many examples can be found in the review paper [31], see also the works |29, [17|, 



18| |16| , [19|j . Here we describe very briefly its main steps. 

A Lax matrix for a given integrable system is a matrix L(u) dependent on a 
"spectral parameter" ugC such that its characteristic polynomial obeys two con- 
ditions 

(i) Poisson involutivity: 

{det(L(u) - v ■ 1), det(L(w) - v ■ 1)} = , Vit, u,v,v G C; 

(ii) det(L(u) — v ■ 1) generates all integrals of motion Hi . 

A Baker- Akhiezer (BA) function is the eigenvector 

L(u) f(u) = v(u) f(u) (3.3) 

of the Lax matrix L(u), provided that a normalisation of the eigenvectors f(u) is 
fixed 

n 

a-f = Yl < u ) /*(«) = 1 » ( /(«) = (A(«). • • • » /»(«))' ) • ( 3 - 4 ) 
i=i 
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The pair (u, v) can be thought of as a point of the spectral curve 

det(L(w) -v ■ 1) = 0. 



(3.5) 



The BA function f(u) is then a meromorphic function on the spectral curve. 
The recipe for finding an SoV is simple: 

The separation variables Uj are poles of the Baker- Akhiezer function, provided it is 
properly normalised. The corresponding eigenvalues Vj of L(uj), or some functions 
of them, serve as the canonically conjugated variables. 

It is easy to see that the pairs (iij, Vj) thus defined satisfy the separation equations 
(|3.2| ) for $j = det(L(uj) — Vj ■ 1). The canonicity of the variables (v,j,Vj) should be 
verified independently. No general recipe is known how to guess the proper (that 
is producing canonical variables) normalisation for the BA function. In many cases 
the simplest standard normalisation, 



a[u) 



a = (0,0,...,0,l) 



(3.6) 



works. In other cases the vector a may depend on the spectral parameter u and the 
dynamical variables (x,p). We shall refer to such normalisation as a dynamical one. 

From the linear problem (|3.3|) and normalisation ( |3.4|) we derive that a ■ L k f = 
v k , k = n — 1, hence, 



/ = f(u) 



( a 
a ■ L{u) 



( 1 \ 

v 



V v 71 - 1 ) 



(3.7) 



Another useful representation of the eigenvector f(u), which can be directly verified, 
is as follows: 



{L{u) - v ■ 1)A 
(S.(L(u)-v-l)*) k ' 



Wk = 1 n . 



(3.J 



where the wedge denotes the classical adjoint matrix (matrix of cofactors). 

To derive equations for the separation variables, let = res fi(u) and Vj = 

). Then from ( p.3|) -( ^4|) we have the overdetermined system of n + 1 linear 



v\u q 



homogeneous equations for n components f\ of the vector f^: 

f L( Uj ) /W = Vj /W , 

The pair (u,v) = (uj,Vj) is thus determined from the condition 



(3.9) 



rank 



L(u) -v ■ 1 



1 . 



(3.10) 
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Finally, the condition fl3.10| ) can be rewritten as the following vector equation: 

a- (L(u) -v- 1) A = 0. (3.11) 

One can eliminate v from ( |3.11|) to get the equation for u/s in the following way. 
From the linear system ( |3.9| ) it follows that a ■ (L{uj)) k f^' — 0, k = 0, . . . , n — 1, 
so that (because f^' is not a zero vector) the following determinant has to vanish 
on the separation variables uf 



B(u) = det 



/ a \ 
a ■ L{u) 

\a-L n -\u) j 



(3.12) 



The formula ( |3 . 1 2| ) for the separation variables appeared already in |24| (see also 
0) in the case of standard normalisation: a = olq (|3.6| ) (see, for instance, formula 
(22) in 123). 

Notice that the fact, that equations ( |3.11|) and fl3.12j ) are the ones for the poles 
of the BA function, is already hinted, respectively, by the formulas ( |3.8| ) and (|3T7|). 
Also, from equations ( |3.11| ) we can get many various formulas for v in the form 



v = A(u) 



(3.13) 



with A{u) being rational functions of the entries of L(u). Let us describe those 
formulas for A(u) explicitly. 

Define the matrices Lfo\ p = 1, . . . , n, with the following entries: 



-(p) 



E- E 



ii=l 



ip-i- 



J il ,j 



Lip-i,j Li p _ 1 i l 



p = 2,3,...,n, (3.14) 



1p— 1 — 1 

and put = L. These matrices satisfy the recursion relation of the form 

L (p) = L (trL^) - (p - 1) . 

Introduce the matrix B(u) by the formula 

/ 



(3.15) 



B{u) 



a- lW(w) L-\u) \ 
a-LW(u) L-\u) 
\a-L^(u) L-\u) 



(3.16) 



\^3-LW{u) L-\u)) 

Then we have the following statement. 
Proposition 3. 



a-(L(u)-vir = ((-v) n -\(-v) 



\n-2 



B{u) 



(3.17) 
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Proof. The characteristic determinant det(L(w) — v-1) has the following represen- 
tation 

det(L(u) -vl) = (-v) n + Y ±— { tr(L^(u)) . (3.18) 

pi j\ 

The adjoint matrix {L(u) — v ■ 1) A is a matricial polynomial in v of the degree n—1, 

n-l 

(L(u) -vl) A = {-v)^ 1 1 + (-v) 71 ' 1 ^ A {j \u) . (3.19) 

3=1 

In order to find the matrices A^> , substitute ( |3.18|) and (|3.19| ) into the definition of 
the adjoint matrix, 

det(L(u) - v ■ 1) • 1 = (L(u) - v • 1) (L(u) - v ■ 1) A , 

and equate coefficients with the degrees of v. In this way we get the following 
recursion relation for the A^^s 

A {j) = f { tr(L (j) ) - L A^-V (3.20) 

with the initial data 

A (0) = 1 ) A (n-1) = i. tr ( L W) L -l( u ) _ (3,21) 

The matrix ^ L" 1 (cf. (CT)) satisfies the same recursion and the same initial 

values which means that 

A^\u) = ±L^ +1 \u)L-\u). 

■ 

From the system of linear homogeneous equations 

a ■ (L(u) -vl) A = ((-v) n ~\ (-v) n - 2 , . . . , 1) • B(u) = (3.22) 
(cf. (|3.17|) ) we derive that 



(-r = }*, Vfc. (3.23) 
{B A {u)) kj 

The formula ( p. 23 ) gives plenty of different representations for the function A(u), all 



of them being compatible on the separation variables since, because of the equality 

B{u) = det(B(u)) , (3.24) 

the matrix B A (uj) has rank 1. 

To validate the choice of normalisation a(u) it remains, first, to make sure that 
the number of u^s is exactly iV (in some degenerate cases one has to supply a 
couple of extra variables to make a complete set) and, second, to verify (somehow) 
the canonicity of brackets between the whole set of separation variables, namely: 
between zeros Uj of B{u) and their conjugated variables Vj = v{uj) = A{uj). To do 
this final calculation one needs information about Poisson brackets between entries 
of the Lax matrix L(u). 
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4. The separation 



We now proceed with applying the general method to the system in question. For the 
Ruijsenaars model the number N of degrees of freedom coincides with the number 
n of particles and, respectively, with the dimension n of the Lax matrix ( J2.ll ) , so 
we can put N = n in the formulas of the above Section. Let us first prove two useful 
Lemmas. 

Lemma 1. Let q G C, xi- G T>, i = 1, . . . , M , j = 1, . . . , N, be arbitrary constants 
such that 

N 

xf^ = x (mod T) Vi . 

3=1 

Then there exist C G C, yj G V, j = 1, . . . , N, such that 

M N N 

P( u ) = Ci II a ( u ~ x f) = C\\a{u- yj) , Vn G C , 

4=1 j=l j=l 

where 

N 

y j = x (mod V) . 

3=1 

The p(u) can be thought of as cx-function version of the iVth degree polynomial (in 
u) which is represented in terms of its zeros yj. 
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Proof. Let Zj G T>, j = 1, . . . , N, be N distinct constants such that Y,j=i z 
x (modr). Consider the elliptic function p(u) of the form 

pW ^^i4[' (41) 

Any elliptic function can be represented through the ratio of products of cr-functions 
depending on its zeros, yj, and its poles, Zj (ci., for instance, 0), i.e. 

p>)=c nto4, (4.2) 



i= i °K U ~ z h 



where J2j=iVj = J2j=i z j = x (modT). The statement follows if we equate right 
hand sides of (OF) and (OI). ■ 



Consider the Lax matrix L(u) for the Ruijsenaars system 



L(u) =J2 h i ®u(xi - Xj + A) Eij . (4.3) 

i,j=l 
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Lemma 2. For any integer p = 1, 2, . . . , n we have the identity 



{L{u)) P + EW -4t~, a{u + 3X) (L(u)y-i =±K C$ - x 3 + pX) E tJ 

j=i aJ \ A ) a \ u ) i,j=i 

(4.4) 

where the scalars C^f do not depend on the spectral parameter u and are given by 
the formula 

r»(p) _ t i \p-i u h IIfc<; a ( x ik ~ x h) a ( x H ~ x iJ v ( A r\ 



cr(xj — Xj + pA) 
cr(xj — Xj + A) 



p-i 

n 

■? 1 'V fe=l 
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(x^ x^ ) o~(xi k x j) 



a(xi - x ik + A) cr(x ife - + A) 



(i) 



and C)a 

This Lemma does actually say that it is possible to arrange for the degree p poly- 
nomial in L(u) (the left hand side of (|4.4|) ) such that w-dependence of its (ij)-entry 
occurs only through the factor $ u (xj — Xj + pX). This fact reflects some hidden 
internal structure of the Lax matrix L(u) and is essential for further proof of the 
separation of variables. Notice also that the usage of the generalised Cauchy identity 
is very important for the proof of the Lemma given below. 

Proof. Iterating the recursion ( ft.!5|) for the matrix L^ p \u), we get the formula 



l (p) = (-iy-\ p - 1)! L p + ^(-l)P-W Jj: tr (L^) L p - j . (4.6) 

3=1 



Noticing that the traces of the L^' matrices are expressed in terms of the integrals 
of motion (cf. ( fTTSP and (PISQ) 

tlLW=jl ^°±+M (47) 

we have that 

In C\f $ u (x, - x, + pX) = 2$ . 

The right hand side being evaluated with the help of the generalised Cauchy identity 
(|2.7p , we arrive to the statement of the Lemma. ■ 

In order to separate variables in the Ruijsenaars system, first of all we have to 
fix the normalisation vector a(u). The crucial observation is that we can use the 
standard normalisation (|3.6| ). Then we have the following "characteristic equations" 
for the separation variables u = Uj and v = Vj (cf. ( |3.11| )) 

(L(u)-vl) A nk = 0, fc = l,...,n. (4.8) 
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The "a-polynomial" B(u) Q3.12| ) has now the form 



B (u) = det 



/ ... 1 

B n \ ■ ■ ■ L nn 

\(L n ... [L n l ) nn J 



(4.9) 



Its zeros, Uj, are the poles of the BA function f(u) and are the separation variables. 
Let us first verify that we have got the right number of the u/s. 

Theorem 1. a -polynomial B(u) ( \4-9j) hasn — 1 zeros Uj G T> and can be represented 
by the formula 



n— 1 



B{u) = C H$ u ( 



-Un 



(4.10) 



where C does not depend on the spectral parameter u and has the form 



C = (-l)"" 1 hi' 1 


1 


1 

u n,n-l 




r (n-l) 


r (n-l) 
u n,n-l 



(4.11) 



Variables Uj obey the restriction 



n-l 



(4.12) 



i=i 



Proof. Using Lemma 2 we can represent B(u) in the form 

= (-1)™" 1 /i™- 1 x 



(4.13) 



®u{x n l + A) 

Ci? $ u (x nl + 2A) 



Cn,n-1 $ u (x n ,ri-i + 2A) 



^u{x n ,n-l + A) 



Cir x) $„(x nl + (n-l)A) 



C (n_1 ]$ fx 



n 



1)A) 



Then, using Lemma 1, we conclude that the u-polynomial B(u) can be rewritten in 
terms of its zeros in the form ( |4.10| ) where C is given by the formula ( |4.11|) and we 
also have the restriction ( 4.1 2|) . ■ 

To avoid discontinuities when discussing the Poisson brackets it is convenient to 
think of Uj S cLS lying on the torus C/T rather then on T>. 

In the sequel we obtain few statements which are valid for a general Lax matrix 
L(u). Let us introduce the following matrices: 



L(u, v) 



L(u)-v-l 



L A (u,v) := (L(u)-v- 1)' 



(4.14) 



We can express the Poisson brackets of L A (u, v) with L(u, v) in terms of the Poisson 
brackets of L(u, v) with L(u, v). The answer is given by the following Lemma. 
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Lemma 3. 

{L A (u,v),L 2 (u,v)} = A^ 1 (tri [ L A (u, v) {Lx(u, v), L 2 (u, v)} ] 

~ L i(u, v) {Li(tt, v), L 2 (u, v)}) L A (u, v) , 

{Li{u,v),L$(u,v)} = A 2 1 (tr 2 [ L 2 (u, v) {L^u, v), L 2 (u, v)} } 

- L 2 (u, v) {Li(«, v), L 2 (u, v)}) L 2 (u, v) , 

where L^u, v) = L(u, v) ® 1, L 2 (u, v) — 1 ® L(£t, ?)), L A (u, u) = L a (m, t> ) ® 1 etc, 
Ai = det(L(u,v)), A 2 = det(L(u,v)) and tri^ means trace in the first, respectively, 
the second space of the tensor product of two spaces and is defined by the rule: 

tri[AiB 2 ] = tr x L4 (8) B] := tr(A) (1 ® S) = tr(A) 5 2 , (4.15) 

tr 2 [AiS 2 ] = tr 2 [A ® 5] := tr(B) (A ® 1) = tr(5) A x . (4.16) 



Proof. The matrix L(w, v) and its classical adjoint L A (u, v) by definition satisfy the 
relation 

L A (u, v) L(u, v) = L{u, v) L A (u, v) = Ai • 1 . (4.17) 
Differentiating this formula with respect to a parameter t and using the formula 

one obtains (cf. (1.45)-(1.47) from 0) 

dL A _ L A tr (l a I L) - L A (i L) L a 



(4.18) 



dt Ax 
From which we have the following derivatives in the component-wise form: 

flfA rA rA rA rA 
dL pq A 1 
Now, using the derivation property of the bracket, 

{Lij(u,v),L H (u,v)} =J2 ^r~T~~\ { L pq(u,v),L M (u,v)} , 

pq OL pq [U,V) 
pq OLi pq [U,V) 

we verify both statements of the Lemma by substitution and straightforward calcu- 
lation. ■ 

From the involutivity of the characteristic polynomials of a Lax matrix L, A x 
and A 2 , we have the equality: 

= {Ax ■ 1 <g> 1, A 2 • 1 ® 1} = {L 1 (u,v)L A (u,v),L 2 (u,v)L A (u,v)} 

= Lt L 2 {L A , L A } + L x {L A , L 2 } L A + L 2 {L u L A } L A + {L h L 2 } L A L A . 



12 



Hence, using Lemma 3, we can get from here an expression for the bracket of L A 
with L A in terms of the brackets of L with L. 

Lemma 4. 

{L$(u,v),L$(u,v)} = A^ 1 A 2 1 ( L A L 2 {L x , L 2 } — tr x [L A L 2 {L x , L 2 }\ 

-tv 2 [L A L A {L^L 2 }}) L A L A . 

Suppose now that a Lax matrix L(u) satisfies the quadratic (dynamical) (r, s)- 
bracket, then we have the following statement. 

Lemma 5. Let a Lax matrix L(u) satisfy the quadratic (r, s)-bracket of the form 
{Li(m), L 2 (u)} = L x L 2 r + - r_ L x L 2 + L x s + L 2 - L 2 s_ L x (4.19) 

where 

r + — r_ + Si — s_ = , Pr ± P = — r+i , V s+V = s=\ 

Here V is the flip in tensor product of two spaces, i.e. V (A® B)V = B <g> A. Then 
the matrix L A (u, v) = (L(u) — v ■ 1) A obeys the bracket of the form 

{L A ,L A } = {r + -ti x r + -ti 2 r + )L A L A -L A L A (r„-tr x r„-tr 2 r-) (4.20) 
+ L A (s + - tri s + - tr 2 s+) L A - L A (s_ - tr x s_ - tr 2 s_) L 2 
+ V Ar 1 [(L A (r + - s _) - tn[L A (r + - *_)]) L A L A 

-L A L A ((r_-s + )L A -tr x [(r„-s + )L A })} 
+ v A, 1 [(L 2 > + + s + ) - tr 2 [L A (r + + s + )}) L A L A 

-L A L A ((r_ + S -)L A - tr 2 [(r_ + S -)L 2 A ])] . 



Proposition 4 ( ||21| , |32|| ). TTie Lax matrix ( 2.11\) of the Ruijsenaars model satisfies 
quadratic (r, s)- algebra where (r, s) -matrices can be chosen to be as follows: 



r + = a — b + c — d . 
where 



r_ = a + d . 



s + = b + d . 



c — d . 



d 



§u-u{xi - Xj) Eij <g> Eji + ((u -u)^2 E kk ® E kk 

i^j k 

®u(xi ~ xj) ® En + C(u)J2 E kk ® E kk , 

i^j k 

$u(xi - Xj) En ® £?y + C(u) E kk ® ^fcfc , 

^ C(x, - Xj) /•;„ (g> £jj . 

i+3 



(4.21) 

(4.22) 
(4.23) 
(4.24) 
(4.25) 
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Notice here that one needs to use three algebraic relations ( |2.5[ )-(0) and ( |2.4| ) for 
the function $ to verify this (r, s)-structure (cf. [[2l|j ). 



Separation variables (u, v) = (uj, Vj), j = 1, . . . , n — 1, for the Ruijsenaars model 
are implicitly defined by the following system of equations 

{L A {u,v)) nk = {L{u)-v -1)^ = 0, k = l,...,n, (4.26) 

where L(u) is the Lax matrix ( |2.11| ). The Poisson brackets for these new variables 
are generally given by the expression: 



{Ui,Uj} {u U Vj} 
{Vi,Uj} {Vi,Vj} 



(M ; w) x (4.27) 



f{(L A (u,v)) nk ,(L A (u,v)) nk } {(L A (u,v)) nk ,(L A (u,v)) nl }\ t a 

X I {(L>, v)) nh (L A (u, v)) nk } {(L>, v)) nh (L A (u, v)) nl } J , ^ kl) 

where it is assumed that k ^ I, the condition means substitution of the form 

(U,v) = (Ui,Vi) 
(U,V) = (Uj,Vj) 

and matrices Ai are defined as follows: 

/ d(L*(u,v)) nk d(L A (u,v)) nk \ 
M m - k l := d ( L 4u,v))„, d(L^u.v))^ • ( 4 - 28 ) 



Aij 



du dv 



(u,v) — (u m} vm) 



Theorem 2. The separation variables (uj, Vj), j = 1, . . . , n — 1, for the Ruijsenaars 
system, defined by the system of equations ( \4-20{ ), possess the following Poisson 
brackets: 

(i) {ui,Uj} = {ui,Vj} = {v h Vj} = 0, i^j, 

(ii) {v j ,u j } = v j . 



Proof. Generically the matrix M. m - k i ( |4.2<j| ) for k ^ I is invertible which means that 
in order to prove the statement (i) we have to show that 

{(L A (u,v)) nk ,(L A (u,v)) nl } lA =0, VM = l,...,n, 

I Ay 

when i j. The latter fact follows from the Lemma 5 when we substitute in the 
right hand side of ( |4.20| ) the (r, s)-matrices from the Proposition 4 and put in both 
sides (w, v) = (ui,Vi), (u,v) = (uj,Vj), i ^ j. Indeed, using the definition fl4.26|) , we 
get then the expression of the form 

{(L A {u, v)) nk , (L A (u, v)) nl } u (i + j) (4.29) 

= J2 {{a-b + c- d) np!nq )\ A .. (L A (ui,Vi)) pk {L A {u h Vj)) q i 

pq 

+ £ ( - ^M)^^^ {a _ b)pqnr \ (L A ( Uj , Vj )) rl 

pqr ' A ij 

2^1 dct(L(u,v)) \ U ^~ c )nr,pq), \^h"i))rk- 



pqr' 
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Each of the above three terms is equal zero, the first one when simply inspecting 
the inputs from the matrices a, b, c, d; the latter two because the simple zero in the 
denominator is cancelled by a double zero in the numerator. 

In order to prove the statement (ii) we take % = j in ( [4.27 ) to get 



{ Uj , vj} det Mj-u = {(L A (u, v)) nk , (L A (u, v))ni}\ MHujiVj) , (4.30) 
where we recall that k ^ I. Hence, we have to show that 

- Vj det MyM = {(L A (u,v)) nk , (L A (u,v)) nl } hu . (4.31) 



To calculate the right hand side of Q4.31|) we use the Proposition 4 and take the limit 
u — > u in the (r, s)-bracket ( |4 . 1 9| ) . Using the derivation property of the bracket and 
substituting u = Uj, v = Vj we then conclude that the only non- vanishing term in 



the right of (|4.31|) has the following form: 



{(L A (u,v)) nk ,(L A (u,v)) nl } lM ^ 



= ff d(L*(u,v)) nk d(L*(u,v)) nl _ d(L A (u,v)) nl d(L*(u,v)) nk \ d(L(u,v)) ps 

3 Z-^i \ \ d(L(u,v)) pr d(L(u,v)) ra d(L(u,v)) pr d(L(u,v)) rs J du 
prs 

On the other hand the determinant of M.j±i can be evaluated making use of its 
definition flOg) and expressing the derivatives by u and v in terms of those by 



(L(u,v)) pg . Then we have the following formula for the left hand side of ( [4.311 ): 



- Vj det Mj. k i (4.32) 

\ - ( ( d{L*{u,v)) nk d{L A (u,v)) nl d(L*(u,v)) nt d(L*(u,v)) nk \ d(L(u,v)) pa ' 



V J Z-^i \\ d(L(u,v)) ps d(L(u,v)) rr d(L(u,v)) pa d(L(u,v)) rr J du , . 



Straightforward calculation, using ( 4.18|) and the fact that the matrix L A (uj,Vj) 



has rank 1 shows that these two expressions are equal to each other (cf. here the 
proof of the analogous Theorem 1.3 from |J establishing the Poisson brackets for 
the separation variables for the sl(n) Gaudin magnet obeying the simplest linear 
r-matrix algebra with the rational r- matrix). ■ 

Theorem 3. The variables {uj,yj := log(u,)) ; j = 1, ... ,n — 1, together with the 
variables (X,P) describing the "motion of the center- of -mass" , 

n 

X:=x n , P := \og(H n ) = ]T Pj , (4.33) 

i=i 

constitute the complete canonical set of new (separation) variables. 

Proof. The bracket {P, X} = 1 is easily seen, so, in addition to the statements of 
the Theorem 2, it is only left to check that 

{P, Uj } = {P, Vj } = 0, j = 1, . . . ,n - 1 , (4.34) 
{X,u j } = {X,v j } = 0, j = 1, . . . , n - 1 . (4.35) 
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The equalities ( f4.34| ) are trivial since (uj,Vj), j = 1, . . . ,n — 1, are denned by the 
equations (L(u) — v ■ l)^ k — 0, k = 1, . . . , n, and entries of the matrix L{u) depend 
only on differences Xi — Xj , therefore 

{P, (L(u)) y } = 0, V/../ 1 n. 

For the brackets in (|4.35|) we have the following expression (k ^ I): 



f{X, (L\u,v)) nk } 
{{X, (L A (u,v)) nl } 



(4.36) 



\(u,v) = {uj,Vj) 

The vector on the right of ( |4.36| ) is equal to zero since Vfc = 1, . . . , n 



{X, (L A {u,v)) nk } { 



I U = Uj 
V = Va 



E 



T A T A J A T A 
Ll nk Ll qV~ L 'np L 'qk 



dct(L{u,v)) Wn P 



(L pq + v5. 



pq) 



0. 



I U = Uj 

V = Va 



The equalities (|4.35|) follow because the matrix M.j-ki is nondegenerate. ■ 

The proved SoV for the A n _i (n-particle) problem with the standard normalisa- 
tion vector ct = (0,0,...,0,1) actually implies another SoV for the A n _ 2 problem 
with the non-standard normalisation vector d^: 



<&u{i — %\ + A), . . . , - x n ^ + A) ) 



(4.37) 



if we choose £ = x n . Let us demonstrate this explicitly. 

Let us take the Lax matrix (|2.11 ) for the n-particle system 



L ( u ) = E ^ § u {xi - xj + A) Eij . 

If we remove the last (nth) row and the last column from this Lax matrix then we 
get the following (n — 1) x {n — 1) matrix 



/ 



LTniu) :- 



hi $ tt (A) 



\/in-l $u(Xn-l ~ Xl + A) 



hi - X n _! + A) 



h n -l $u(A) 



(4.38) 



which is the Lax matrix for the integrable system with n — 1 particles with the 
Hamiltonian 



H[* n) = $ U (A) y £h i = $ U (A) £ e Pl n 



n-1 



n-1 



i=l i=l 

Under the simple canonical transformation, 



e Pl -> e Pl 



Cr(Xj - x n - A) 



,Tv ► Xi 



% = h 



a(xj - x k - A) 
a(xi - x k ) 



(4.39) 



, n — 1 



(4.40) 
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the system fl4.39p turns into Ruijsenaars' system with n — 1 particles. This 1-degree- 



of-freedom-less system obviously inherits the non-standard SoV with the dynamical 



normalisation ( |4.37| ) from the standard one (with do) for the system with n degrees 
of freedom. Indeed, to see this, it is sufficient to note that the separation variables 
(uj,Vj), j = 1, ... ,n — 1, for both systems are defined from the intersection of two 
spectral curves: 

f det(L(u)-i;.l) = 0, , ■ 

\ det(L*Z(u) - v • 1) = . 1 ' 

In other words, the condition of the standard SoV for the first problem, 

rankf . "° \ = n - 1 , (4.42) 



L{u) -vl. 

implies the following condition of SoV for the second problem: 



where d\{u) is given by ( [4.37] ). 



Procedure shown above, on how to connect the standard normalisation vector do 
and the alternative one, di, does obviously reflect an embedding, gl(n — 1) C gl(n), 
of one problem into the other. In other words (and it is true in general, for any 
integrable system of A n type), one always has a free choice, namely: to include or 
not to include the "center-of-mass variable", X, and its conjugate one, P, in the 
complete set of separation variables. 



5. Generating functions 



In this Section we derive the explicit formulas for SoV in the simplest cases: n = 3 
with the standard normalisation ( |3.6| ) of a, and n = 2 with the dynamical normali- 
sation (|4.37|) (we skip the trivial case of the purely coordinate SoV x\p — *• X\ ± x 2 
for the 2-particle problem). Since the both cases are treated in very much the same 
manner as their trigonometric prototypes, see, respectively |31| and | 19[| , we present 
only the main formulas here, omitting the details of the calculations. 

Let us start with the n = 3 case. Following |I8| define two functions Ai(u) and 
A 2 (u) by the formulas 



(L(u) - A k 



'3,3-fc 







k = 1,2. 



or explicitly, 



AM 



Lkk — 



L-tk L 



3k J^k,3-k 



e Pk a k (u) 



a k [u) 



a(u + 2X + x 3 - z 3 _ fc ) a(x k - x 3 - k - A) 
a(X) a(u + A + x 3 - x 3 _ k ) a(x k - x 3 _ k + A) 



k = 1,2. 



k = 1,2 



(5.1) 

(5.2) 
(5.3) 
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The separated variables Uj are defined from the equation 

A 1 (u J ) = A 2 { Uj ) (5.4) 

which is equivalent to the equation B(uj) = since 

B(u) = hj $„(x 3 — Xi + X) $ u (x 3 - x 2 + A) tr(A) (A 2 - A x ) 

and has two roots u± >2 G P. From the easily verified invariance of the ratio 
a\{u) / a 2 {u) under the transformation u t— > x\ + x 2 — 2x 3 — 3A — u it follows that 

U\ + u 2 = xi + x 2 — 2x 3 — 3A (modT) , (5.5) 

which agrees with ( 4.12 ). The conjugated variables Vj = e w are defined as 

Vj = Aiiuj) = A 2 {uj) (5.6) 

or, equivalently, through four equations 

Vj = e Pk a k (uj) , j, k G {1, 2} , (5.7) 

for four variables Mi, u 2 , v i, t> 2 . By virtue of the Theorem 3 the variables (ui, u 2l X; 
Vii Vit P) are canonical. The generating function of the separating canonical trans- 
formation M is most conveniently expressed in terms of another set of canonical 
variables 

x + = x\ + x 2 — 2x3 , X- = x\ — x 2 , X = X3 , (5.8) 

P±=i(Pi±Pa)> P = Pi+P2+P3, (5.9) 
u± = u 1 ±u 2 , y± = \(yi ± y 2 ) . (5.10) 
We shall need a a-generalisation of the Euler dilogarithm function, 

U 2 {z)= riog(sin(C))rfC, 
J o 

which we define as 

S(*):= / tog(a(C))dC. (5.11) 

Notice that this function was introduced in and has been used to construct 
the Lagrangian function of the integrable map which is a time-discretisation of the 
Ruijsenaars system. Using the product expansion for the Weierstrass sigma-function 
(g = exp(*7r^)): 

, , 2ui x /itz\ ~ /l-2g 2 ™cos(5f)+g 4 -\ 



cf. ||33|| , we can express the function S in terms of the following function 

oo z k 

Li 3 (*;g) := g (1 _ gfc) k2 > l?l <!> kl <!■ ( 5 - 13 ) 



fc=l 
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Notice that similar, but different from ( |5.13| ), g-deformations of the Euler (di-) 
trilogarithm have been proposed in the review article [11]. In terms of ( 5.13|) we 
obtain 

rjiZ 3 ( /2aV x ~ 2k 



™ = ik + W?) + %lTw)' (5 ' 14) 

+ i ^(u 3 ( q h^)-U 3 ( q H- l ;e 



IT 

where t = exp(iriz/uji). This series representation converges for \q\ 2 < \t\ < \q\~ 2 - 
Let 

C{y\ x, y) := S(u + x + y) + S(v — x + y) + S(v + x — y) + S{u — x — y) . (5.15) 

The generating function F(y + , X-\. J U— j X— ) of the canonical transformation from 
(x±,p±) to (u±,y±), satisfying the defining relations 

OF OF OF OF 

a — = P+> = U +' q — = P-, o — = -2/-» ( 5 - 16 ) 

ax + ay + ax_ aw_ 

is given then by the expression 

(A X— u \ 
2 ; T't) (5 ' 17) 

+ 5(A-x_) + 5(A + s_). 

The case n = 2 with the normalisation 

«i = ( $„(£ - xi + A), $ u (£ - x 2 + A) ) 

(cf. ( |4.37| )) is treated similarly to its trigonometric prototype []19 . Having intro- 
duced the functions A\{u) and A 2 (u) by the formulas 

(a r (L(n)-4) A ) 3 _fc = 0, fc = l,2, (5.18) 

or explicitly, 

A i = L " - fcfeS} L ^ = e» fll (tt) , (5.19) 

A * = ^ - i:tT + x) L * = eP2 > ( 5 - 2 °) 

„ /.A _ v(u+£,+2\-x a _ k )<T((,-x k )<r(x k -x z _ k -X) JL _ i 9 / p: o 1 \ 

- (T (A) ( 7(«+C+A-a:3- fe )a«+A-x fc ) CT (a ;fc - a; 3_ fc +A) ' ^5^5 ^.ZXj 

one proceeds as above with the only difference that the relation ( |5.5| ) is replaced by 

u x + u 2 = xx + x 2 - 3A - 2£ (mod T) (5.22) 

and the variables x± are defined now as x± = Xi ± x 2 - The resulting expression for 
F(y+,x+;u-,xJ) is 

F = y + (x + -3\-2£) + x + logo-(A) (5.23) 

-r(~ — —\-r(- x +~ x c x -\ 

\2' 2 ' 2 J 1,2' 2 ^' 2 J 
+ S(A-x_) + 5(A + x_). 
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6. Nonrelativistic limit to the Calogero-Moser 
system 



The nonrelativistic limit is obtained by letting A — > while rescaling the mo- 
menta Pj := iXpj/g, jet, and making the canonical transformation pj := pj — 
ig J2k=£j C( x j ~ x k) such that hj — > 1 + iXpj/g + 0(A 2 ) in (|2.11|) . The (r, s)-matrix 
structure is linear in that limit since the L-matrix behaves as 



L(u) - (A" 1 + ««))• ! + 0(A) 



(6.1) 
(6.2) 



£ ( u ) : = &i E n ®u(xj - x k ) E jk . 

3 3+k 

The ^-matrix ( |6.2| ) is Krichever's [I2~ Lax operator for the elliptic Calogero-Moser 
system with the Hamiltonian 



(6.3) 



3=1 



Proposition 5 ( [|30|| ). The Lax matrix £{u) ( \6. 3j ) of the elliptic Calogero-Moser 
system satisfies linear (r, s) -algebra of the form 



{£i(u),£ 2 (u)}=[£ 1 ,r] + [£ 2 ,s} 



where 



r = a + c . 



a — b . 



-VrV\ 



(6.4) 
(6.5) 



(see ^j.22j , ^.2^) , ^2^( j), and [. , .] means matrix commutator. 

The SoV for the elliptic Calogero-Moser system follows, in principle, by taking 
limit A — > in the corresponding formulas describing SoV for the Ruijsenaars system. 
Although, because this limit is not so simple and straightforward, we prefer to do it 
independently, repeating the steps for proving main statements for the Ruijsenaars 
system in Section 4. 

The normalisation vector is the same: 



a(u) 



a o = (0,0,...,0,l) 



We have now the following characteristic equations for the separation variables u 
Uj and v = Vj 

0, i 



('(«) - « " l)n* 

The zeros of the er-polynomial b(u) 



( 



b{u) := det 



\{£ 





n-l> 



Ml 



K = 1, . . . , n . 

' / ) 

""nn 
{£ )nn ' 



(6.6) 



(6.7) 
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give us separation variables Uj. 

Theorem 4. a -polynomial b(u) ( \6. 1\ ) has n—1 zeros Uj G T> and can be represented 
by the formula 

n-l 

b(u) = C 1] ®u{-Uj) (6.8) 

3=1 

where C does not depend on the spectral parameter u. Variables Uj obey the restric- 
tion 

n—1 n—1 

Yl U J = Y ( x i - x n) (mod r) • (6.9) 

3=1 3=1 



Proof. From the limit ( j6.1|) and the definitions of B{u) and b(u) we conclude that 

B(u) = b{u) + 0(X) . 

Both B{u) and b(u) are cr-polynomials in u and, since the degree of such a polynomial 
must not change with the analytical continuation of the parameter A, b(u) has the 
same degree as B(u) does. Moreover, now the separation variables have to obey the 
restriction ( |6.9|) , the one being the limit of the corresponding relation ( 4.12[) . ■ 

Let us introduce the following notations: 

£(u, v) := £(u) - v ■ 1 , t{u,v) := (£(u) -v 1) A , (6.10) 

and also Ai = det(£(u, v)) and A 2 = det(£(u, v)). Suppose now that a Lax matrix 
£(u) satisfies the linear (dynamical) (r, s)-bracket (|6.4j) , then we have the following 
statement. 

Lemma 6. Let a Lax matrix £(u) satisfy the linear (r, s)-bracket of the form 

{£ x {u),£ 2 {u)} = [£ x ,r] + [£ 2 ,s] , s = -V r P { _ a . (6.11) 

Then the matrix £ A (u, v) = (£(u) — v ■ 1) A obeys the bracket of the form 

= Ar 1 [(^s-tri[^a])^^-^^«-tn[ a ^])] (6.12) 
+ A 2 1 [ ( t 2 r - tr 2 [ t 2 r] ) t x t 2 - t x t 2 ( r t 2 - tr 2 [ r t 2 ] ) ] . 



Theorem 5. The separation variables (uj,Vj), j = 1, . . . ,n — 1, for the elliptic 
Calogero-Moser system, defined by the system of equations ( |6'. 6 ), possess the follow- 
ing Poisson brackets: 



(l) {Ui, Uj} = {Ui, Vj} = {Vi, Vj} = , 

(ii) {v j ,u j } = l. 
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Proof. In analogy with the proof of the Theorem 2 we have to show first that 

{(£ A (u, v)) nk , (£ A (u, v)) nl } Uij = , Vfc, I = 1, . . . , n , 
when % ^ j. We have from Lemma 6 and Proposition 5 that 

{{£ A (u,v)) nk ,{£ A (u,v)) nl } { 



{i^j) (6.13) 



+ 



[Af 1 (^s-tr 1 [^s])^^ + A 2 1 {£ A r-tT 2 [£ A r])£ A f 2 

\ - \ {^(u,v)) np ^(u,v)) qk -{^{u,v)) nk {l^{u,v)) qp , ,N 1 (pA( XX 

pgr 1 ij 

|" (l A (n,£ 



pqr 



))np(<? A --(l A (t^)) ni (l A 

det(^(n,C)) 



{t{u u v t )) 



rk • 



These two terms in the right hand side have the same form as latter two in ( [4.29| ) 
and, again, they are equal to zero since in both expressions the simple zero in the 
denominator is cancelled by a double zero in the numerator. 

The matrix of derivatives Ai instead of ( |4.28| ) has now the form 



/ d(i A (u,v)) nk d(£ A (u,v)) nk 
KA , , — I du dv 

\ du dv 



(6.14) 



(u,v)= (um,Vm) 



In order to prove the statement (ii) we have to show that 

- det Mj. kl = {(£ A (u,v)) nk , (i A (u,v)) nl } l{u , 



(6.15) 



where k ^ I. Again, the right hand side of ( p,15| ) can be evaluated by first taking the 
limit u — > u in the (r, s)-bracket of the Proposition 5 and then using the derivation 
property of the bracket. We derive the following expression 



{(£ A (u,v)) nk ,(£ A (u,v)) nl } l(v 



z^\\ a(£(u,v)) pr 



A '« k d(t A (u,v)) nl _ d(t A (u,v)) nl 8(t A (u,v)) nk \ d{l{u,v)) pl 



prs 



d{£(u,v)) r 



d(e(u,v)) pr d(£(u,v)) r 



du 



On the other hand the determinant of M.j- k i (cf. (|4.32| )) has the form 
- det Mj-ki 



(6.16) 



E 

prs 



d(e A (u,v)) nk d(e A (u,v)) nl d(e A (u,v)u d(e A (u,v)) nk \ d(e(u,v)) ps 



d{t{u,v)) p3 d(£(u,v))rr d(l(u,v)) ps d(e(u,v)) rr 



du 



(u,v) = (u i ,v i ) 



Such two expressions are equal to each other by the reasons pointed out in the end 
of the proof of Theorem 2. ■ 
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Theorem 6. The variables (uj, Vj = yj), j = 1, . . . , n — 1, together with the variables 
(X, P) describing the motion of the center- of -mass, 



X := x n , P := tr i(u) = ^Pj 



(6.17) 



constitute the complete canonical set of new (separation) variables. 

Proof repeats the proof of the Theorem 3. ■ 

Consider now nonrelativistic limit of the generating functions F (|5.17 ) and fl5.23Q 
in the two simplest cases. In analogy with calculations in the previous Section, for 
the case n = 3 let us define two functions A\(u) and ^(w) by the formulas (|5.1|) , or 
explicitly, 



AM 



Lkk ; = Pk + W a k {u) 



k = 1,2 



(6.18) 
1,2. 

(6.19) 

The generating function F(y + ,x + ;u-,X-) (cf. formula (7.12) in PH ) is then given 
by the expression 



-^3,3-fc 

Ofc(w) = C(«) + C( x 3 ~ x k) + C( x k - ars-fc) - C(« + x 3 ~ x s-k) , k 
The ±-variables are defined by (|5.8| )- (|5.10|) and we have the restriction 

u + = x + (rnodT) . 



F = y + x + + ig log 



a 



a 



(6.20) 



Similarly, in the case n = 2, the normalisation vector is taken as follows: 

ai = -a*), $ u (£-x 2 )). 

Introduce the functions Ai(u) and A 2 (w) by the formulas ( |5.18| ), or explicitly, 

A i = L n " fcfer = Pi + ig aM , (6.21) 



2 = £a-£^£ 21 =J> 2 +iWtt) 



(6.22) 



afc(«) = ((«) + ((£ - ^fc) + C( x k - x 3 - k ) - ((u + £ - x 3 _ fc ) , fc = 1, 2 . 
The variables x± are defined in this case as x± = x% ±x 2 and we have the restriction 

u+=.x+-2£ (rnodT). (6.23) 

The generating function F(y + • 3j _i_ . U - 3j ) has the following form 



F = y + x + + ig log 



a: + +x_-2g 



(T 



-2€ 



(6.24) 
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7. Concluding remarks 



We have performed the separation of variables for the classical n-particle Ruijsenaars 
system. If we replace the a-function u{x) in all the above formulas by sin(a;) (sinh(x)) 
or by the identity function: x —>■ x, then we get all the above statements valid for 
the cases of trigonometric (hyperbolic) or rational Ruijsenaars system, respectively. 

We have found the explicit generating function F(u\x) of the separating canonical 
transform in the cases of two and three particles. It is a challenging problem to 
obtain such a function for n > 3 in any explicit form. What is also a problem 
for possible further studies of this integrable system is to produce a quantum SoV, 
i.e. to find the corresponding kernel M.%{u\x) of the quantum separating integral 
operator and related integral representation for eigenfunctions of the quantum 
integrals of motion Hj (cf. © O. S H> 13 ) • 
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